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We implement the concept of complete local symmetry in lossy acoustic waveguides. Despite the 
presence of losses, the existence of a spatially invariant current is shown theoretically and observed 
experimentally. We demonstrate how this invariant current leads to the generalization of the Bloch 
and parity theorems for lossy systems defining a mapping of the pressure field between symmetry 
related spatial domains. Using experimental data we verify this mapping with remarkable accuracy. 
For the performed experiment we employ a construction technique based on local symmetries which 
allows the design of setups with prescribed perfect transmission resonances in the lossless case. Our 
results reveal the fundamental role of symmetries in restricted spatial domains and clearly indicate 
that completely locally symmetric devices constitute a promising class of setups, regarding the 
manipulation of wave propagation. 


I. INTRODUCTION 

Symmetries play a fundamental role in our search 
for the physical laws in nature. Significant theoretical 
and phenomenological advances have been obtained us¬ 
ing models characterized by global symmetries, including 
the existence of forbidden bands of wave propagation in 
infinite periodic structures (Bloch theorem) [T] and the 
existence of even or odd eigenstates in problems that in¬ 
volve potential with global inversion symmetry (parity 
theorem) [2]. However, global symmetry and its conse¬ 
quences constitutes an idealized scenario, requiring often 
a sophisticated preparation procedure to be experimen¬ 
tally probed. Opposite to this, a frequently encountered 
scenario in nature concerns complex systems with a vari¬ 
ety of symmetries valid in limited spatial domains of the 
entire embedding space [3]. For their characterization 
we apply here the term local symmetries and their emer¬ 
gence can be attributed to mechanisms of global symme¬ 
try breaking leading to new symmetries (remnants) at 
different scales. 

Indeed, such spatially localized symmetries can 
be intrinsic in complex physical systems like large 
molecules HE], quasicrystals EHHj, or even in partially 
disordered matter [9]. Furthermore, since technological 
advances often require tailored structures suitable for cor¬ 
responding applications, local symmetries can be present 
by design in multilayered photonic devices [TOHT^ , semi¬ 
conductor superlattices [13] and magnonic systems [14]. 
Acoustic [EHn] and phononic [T31 - BU| structures have at¬ 
tracted increasing interest in this direction. Despite the 
fact that systems belonging to the aforementioned classes 
have been extensively investigated pTH35] , the impact of 
their local symmetries has not been explored thoroughly. 

The natural step beyond the above ascertainment 
would be to investigate whether -and how- these lo¬ 


cally existing symmetries carry imprints of the equiva¬ 
lent global symmetries and in turn, by which procedure 
these could be used to derive structural characteristics of 
the scattered wave field. Recently, we developed a rig¬ 
orous formalism [26] in order to describe the impact of 
local symmetries in one-dimensional, lossless scattering 
setups. To this aim we employed a generic wave mechan¬ 
ical framework -equally applicable to acoustics, photon¬ 
ics and quantum mechanics- described by the generalized 
Helmholtz equation A''{x) + U{x)A{x) = 0, where A{x) 
is the complex wave field and U (x) is the potential which 
describes the inhomogeneities of the medium. For the 
case of acoustic waves, on which this work focuses, A{x) 
corresponds to the pressure field p{x). 

In one dimension, the basic symmetry transformations 
are reflection through a point a and translation by a 
length L, both described by the form F{x) = x = ax + p, 
allowing for a unified description as: (i) Reflection: a = 
— 1, p = 2a, (ii) Translation: cr = 1, p = L. Let us 
assume that U{x) obeys the symmetry U{x) = U{F{x)) 
within an arbitrary domain I? C K. If I? = K then 
the symmetry is global. In any other case the symmetry 
will be regarded as loeal and can be of different kinds, as 
shown in Fig. (for a detailed description of the different 
types of local symmetries see [25] )■ The symmetry of the 
potential leads to a pair of invariant, symmetry-induced 
nonlocal currents, 

Q = ^[aA{x)A'(x) - A{x)A!{x)], (1) 

Q = —[crA*{x)A'{x) — A{x)A!*{x)], (2) 

2i 

which remain constant within F). The term “nonlocal” 
stems from the fact that Eqs. 0, (§ involve both sym¬ 
metry related points x and x. For a detailed derivation 
of Q and Q see Ref. [25]. 
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Especially the current Q has a twofold fundamental 
functionality: (i) it provides a systematic pathway to¬ 
wards the description of discrete symmetry breaking and 
(ii) it generalizes [55] the well-known Bloch [T] and par¬ 
ity m theorems for systems with broken translation and 
reflection symmetry, respectively. 

The symmetry of the potential, U{x) V x G K, im¬ 
plies the commutation of the Helmholtz operator H — 
(P/dx^ + U{x) with the considered symmetry operator 
Op. Interestingly, by setting Q = 0 it is straightforward 
to show that A{x) = OfA{x) = cA{x), which implies 
that the wave field A{x) becomes an eigenfunction of 
the respective symmetry operator Op- This case corre¬ 
sponds to a globally symmetric system, as shown in Fig.[^ 
(a). On the other hand. Fig. (b) illustrates a poten¬ 
tial which is symmetric {U{x) = U{x), V x G I? where 
V is the domain which delimits the potential), though 
the symmetry of the system is globally broken due to the 
asymmetric asymptotic conditions. In this typical scat¬ 
tering situation, even though ['H,Op] = 0, the asymp¬ 
totic conditions prevent the wave field ^(x) from being 
an eigenfunction of Op. Then while Q is spatially con¬ 
stant within the extent of the setup, its value is different 
from zero. Breaking further the symmetry, renders the 
potential nonsymmetric itself C/(x) ^ U(x), Vx G 25, and 
therefore \H,Op] ^ 0. If, however, local symmetries are 
still retained in restricted spatial domains, as 25i, 2?2 in 
Fig. n(c), then Qi, Q 2 are constant within 2?i, V 2 and 
Qi tQ 2 , in general. 

Hence, a vanishing Q value characterizes a system with 
a global, unbroken symmetry, whereas a nonzero Q value 
signifies the breaking of a symmetry. Note, that in the 
extreme case of a total symmetry breaking (e.g. strong 
disorder) Q = Q{x) becomes a function of x, lacking the 
property of constancy, as a result of the lack of any sym¬ 
metric region residing in the potential landscape. In this 
sense, a nonvanishing Q value points to the globally bro¬ 
ken symmetry, whereas its constancy, evinces the pres¬ 
ence of local order and is interpreted as a remnant of the 
corresponding broken global symmetry. This procedure 
allows to track -in a quantifiable manner- the way from 
perfect order and global symmetry to aperiodic systems 
with local symmetries and to disordered systems where 
no local order is present. 

The power of Q as a symmetry breaking measure, was 
confirmed in a recent work m where, having the role 
of the natural order parameter of the system, it allowed 
for the extension of the phase diagram of the well-known 
scattering {VT) symmetric dimer |55|. Further, studies 
of scattering systems involving discrete symmetries re¬ 
vealed intriguing properties [29H3T] of this pathway to 
symmetry breaking. 

In the present work we firstly experimentally verify the 
spatial invariance of the current Q within domains of re¬ 
flection symmetry in aperiodic, completely locally sym¬ 
metric (CLS) [55] acoustic waveguides. The latter can 
be completely decomposed into domainwise locally sym- 
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FIG. 1. (Color online) (a) Globally symmetric system. The 
corresponding symmetry is obeyed V x G K and the wave 
field is an eigenfunction of the symmetry operator, (b) Even 
though the potential itself is symmetric V x G 25, the asym¬ 
metric asymptotic conditions break the global symmetry of 
the system, (c) The potential is not symmetric V x G 25. Nev¬ 
ertheless, symmetry is retained locally in the domains 25i, 252. 
(d) Completely locally symmetric (CLS) setup, consisting ex¬ 
clusively of potential domains with local symmetries. 


metric parts (see Fig. [^(d)). We extend the developed 
theoretical framework in order to incorporate losses and 
find the surprising fact that Q remains invariant within 
the symmetry domains in spite of the presence of losses. 
The parity and Bloch theorems are here generalized for 
inhomogeneous scattering media with losses. As a result 
a mapping of the pressure field values between two sym¬ 
metry related domains is obtained. Focusing on the case 
of parity, we experimentally confirm both the existence 
of the constant current Q as well as the validity of the 
associated mapping. Our aperiodic acoustic waveguide 
is designed according to the construction principle pro¬ 
posed in Refs. [55] [50] for CLS setups with prescribed 
perfect transmission resonances (PTR). The key ingre¬ 
dient allowing for the direct experimental observation of 
the current Q is the ability to measure both the magni¬ 
tude as well as the phase of the pressure field. Therefore 
CLS acoustic waveguides (and corresponding phononic 
structures) represent a useful test bed for newly devel¬ 
oped wave mechanical concepts. 

The paper is organized as follows. In Sec. [^ we de¬ 
scribe the experimental setup and briefly report on the 
theory of wave propagation in acoustic waveguides. In 
Sec.|TH]we describe the procedure followed for the design 
of the device and the way that losses affect the expected 
outcome. Subsequently the existence of the symmetry- 
induced, spatially invariant current Q is proven in the 
presence of losses and we provide the corresponding map¬ 
ping relations which constitute the generalization of the 
parity and Bloch theorems. Sec. jlVj contains the results 
of our analysis, including the comparison of the trans¬ 
mission diagrams between the lossless theory and the ex- 
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periment, the demonstration of the pressure field spatial 
profiles, the experimental confirmation of the invariance 
of Q and the resulting pressure field mapping between 
symmetry related spatial domains. Finally, Sec. |V] pro¬ 
vides our conclusions. 


II. ACOUSTIC WAVEGUIDE WITH 
SIDE-LOADED POINT SCATTERERS 

A. Experimental setup 

In Fig. (a), we show a schematic of the aperiodic 
acoustic structure that was used in this work to experi¬ 
mentally characterize the transmission coefficient and the 
pressure fields of different frequencies. It consists of a 
cylindrical tube with section S = 7ri?^, perforated at the 
positions Xn with holes of diameter dn = 2r„ and length 
£, the thickness of the tube’s wall. The sound source is 
a piezo-electric buzzer embedded in the impedance sen¬ 
sor (see [22] for details) which is placed at one end of 
the main tube. This sensor is used for the calculation 
of the transmission diagram and for the measurements 
of the pressure fields inside the lattice at prescribed fre¬ 
quencies. At the other end of the tube, a handmade 
termination is placed. This termination is a structure 
composed of a purely resistive metal tissue coupled with 
a cavity of tunable length (see appendix of [33] for more 
details). By calculating the reflection coefficient of a 
tube, we found a reflection coefficient less than 5% for 
the frequency range of 100 — 1000 Hz. Thus, we con¬ 
sider this structure as an anechoic termination for our 
further analysis. A 1/2 inch microphone (B&K 4136), 
carefully calibrated, is used, placed after the last hole, for 
the calculation of the transmission coefficient and a 1/4 
inch microphone (GRAS Type 26AC) within the tube for 
field measurements. The frequency range of the applied 
signal is below the first cutoff frequency of the waveg¬ 
uide, fcoi = 4061 Hz, and thus the propagation can be 
considered one-dimensional. 


B. Wave Propagation 

We consider the propagation of sound waves within 
an air-filled cylindrical tube of section S perforated with 
cylindrical holes of section s„ at positions Xn- If the 
geometric characteristics of the holes are much smaller 
than the wavelength of the propagating sound wave, i.e. 
krn 1, the holes can be considered as point scatterers. 
In a section of the tube between two consecutive point 
scatterers, the acoustic wave described by the pressure 
p(x, t) is the solution of the wave equation: 

d‘^p{x,t) _ J_ d'^p{x,t) ^ 

9x2 c2 > 

where c is the speed of sound in air. At each point scat- 
terer (hole), the boundary conditions require the conser¬ 


vation of acoustic flow and continuity of acoustic pres¬ 
sure: 


=-^Vh{Xn,t) (4) 

(5) 

where s„ is the cross section of the hole, v{x, t) 
and Vh{xn,t) are the acoustic velocities inside the tube 
and inside the hole respectively. Taking the time 
derivative of Eq. Q and using the \D Euler equation, 
we obtain 


dp{x,t) 

dx 


dpjx, t) 
dx 


Xn 


pSndVh{Xn,t) , . 

s dt ^ ^ 


where p is the air density. Eqs. and lead to the 
inhomogeneous wave equation | 33 |: 


d'^p{x,t) 
9a; 2 


1 d'^p{x,t) 
c2 dt^ 




( 7 ) 


Eor plane waves the acoustic pressure p{x, t) and the ve¬ 
locity Vh(x,t) are given by: 

p{x,t) = p{x)e~''‘^* and Vh{x,t) = Vh{x)e~'‘^^, 


where uj = ck is the angular frequency. The amplitudes 
p{x) and Vh{x) at x = Xn are related by the impedance 
relation Zn = p{xn)/{snVh[xn))- At each connection be¬ 
tween the waveguide and the hole, the wave impedance 
(which is modified by a change of cross sectional shape), 
the acoustic velocity and therefore the derivative of the 
pressure dp/dx, are discontinuous functions. Using Eq. 
0. the pressure p{x) is described by the following inho¬ 
mogeneous Helmholtz equation: 

+ k'^p{x) = U{x)p{x), (8) 


where 


U{x) ='^S{x - Xn)Un, (9) 

n 

and 

Un = -ik^. ( 10 ) 

= pc/S is the characteristic impedance of the waveg¬ 
uide, Zn is the acoustic impedance of the n*^ hole at 
the point x = Xn seen from the guide. In general the 
hole impedance Zn = Rn — iXn is a complex function 
of the wavenumber, while the real part accounts for 
dissipation. 
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FIG. 2. (Color online) (a),(b) Pictures of the experimental setup, (c) Schematic of the experimental setup, (d) 5-barrier 
correspondence of the acoustic structure in the low frequency regime. Ui,U 2 correspond to the strength of the 5-barrier potential 
and Di, D 2 to the domains of the reflection symmetry decomposition of the waveguide.. 


C. Description of the losses 

The sound waves that propagate in the proposed 
acoustic structure are subject to two different kinds of 
losses. Propagation losses due to viscothermal effects on 
the wall, and point losses due to the scatterers. Vis¬ 
cothermal losses are taken into account by considering a 
complex expression for the wave number. In our case, we 
used the model of losses from Ref. [3S], namely we re¬ 
place the wave number and the characteristic impedance 
by the following expressions 

fc=5^(l + ^(l + (7-l)/x)) (11) 

^c = f(l + ^(l-(7-l)/x)) (12) 

by setting s = r/6 where 6 = is the viscous bound¬ 
ary layer thickness, /r being the viscosity of air, x = y/Pr 
with Pr the Prandtl number, /3 = (1 — i)ly/2, 7 the heat 
capacity ratio of air and r the radius of the considered 
tube. 

Scatterer losses are taking into account in the real part 


of the acoustic impedance of the holes. They come from 
viscothermal losses at the boundaries of the holes wall 
and radiation losses to the outer environment. Thus, the 
impedance of the hole can be expressed as: 

Zn = Rn - iZh tan(fc(£ -|- -|- fo.n)), (13) 

where i?„ = for A:r„ < 1 [36], Zh = pc/s„. 

and to n are length corrections due to the radiation 
inside the principal waveguide m and to the outer en¬ 
vironment [SH] respectively (see appendix for detailed 
calculations). 

Note that for the frequency range of our interest, k{£ + 
ti,n + to,n) "C 1. Thus, ignoring the losses (i?„ = 0) 
and after Taylor expansion of the tangential function, 
Zn ~ —iZh X k{t + ti^n+to,n)- In that case, the potential 
U {x) becomes a series of 6 barriers with strength: 


S{1 + ti^n + to,n) 


(14) 
















































5 


III. ACOUSTIC WAVEGUIDE WITH 
COMPLETE LOCAL SYMMETRY 

A. Construction rules for perfect transmission 
resonances 

The transmission spectra and the control of their prop¬ 
erties in aperiodic, inhomogeneous media is a field of in¬ 
tense study among a large class of wave scattering sys¬ 
tems [5UMT] . The design of our experimental setup has 
been based on the (lossless) construction rules [5SJ|3D] for 
perfect transmission resonances, which provide the tech¬ 
nique for the design of setups which are transparent at 
preselected frequencies. 

In general, isolated perfect transmission resonances are 
expected to occur in reflection symmetric lossless sys¬ 
tems. On the contrary, in aperiodic, nonsymmetric sys¬ 
tems, the transmission spectrum is usually characterized 
by T = < 1, where t is the transmission coeffi¬ 

cient. Nevertheless, the possible complete decomposi¬ 
tions of an aperiodic CLS setup into locally reflection 
symmetric potential parts, has a major impact on its 
transmission properties and particularly on the emer¬ 
gence of perfect transmission resonances. The correspon¬ 
dence between perfect transmission and local symmetries 
has been demonstrated for quantum mechanical |29j and 
photonic [30] systems. Recent results iniiiaiisi, report¬ 
ing on perfect transmitting resonances in globally asym¬ 
metric setups, are unambiguously explained within the 
framework of local symmetries. 

Perfect transmission resonances can be distinguished 
into the following two classes according to the spatial 
profile of \p{x)\ along the device [ 50 ] : 

s-PTRs: Symmetric PTRs are in one-to-one corre¬ 
spondence with the multitude of the different possible 
decompositions of the setup into reflection symmetric po¬ 
tential parts [29]. In this case, the spatial profile of the 
wave field follows the local reflection symmetries of the 
selected decomposition. In order to construct N s-PTRs 
at N preselected frequencies /i, i = N differ¬ 

ent such decompositions should be identihed. Following 
the construction principle developed in [SS] [30] , we can 
find the relevant (5-barrier parameters (strengths Un and 
distances Ln) which allow for the specific perfect trans¬ 
mission properties. Subsequently, these values are used 
to determine the geometric characteristics of the setup, 
e.g. the radius of the holes which play the role of the 
point scatterers, and their distances. 

a-PTRs: The procedure to construct an asymmetric 
PTR requires to divide the setup into two arbitrary (not 
necessarily reflection symmetric) parts TZi and 7^2 ■ The 
first step is to compute the transmission spectra for TZi 
and 7^2 independently and hnd the frequencies of their 
intersections. At these frequencies, a-PTRs are possible 
to occur by hnding the suitable distance between TZi and 
7 ?. 2 , such that the relevant conditions shown in Ref. m 
are satisfied. In this case, the spatial profile of the wave 
field does not follow any (local) symmetry within the 


corresponding setup. 


B. Design based on PTR construction rules and 
minimization of the losses 


Based on the aforementioned construction technique 
we theoretically design a setup which -in the lossless 
case- manifests both kinds of PTRs. Subsequently, the 
emerging geometric characteristics ((5-barrier strengths 
and their distances) are implemented in the fabrication 
of the acoustic waveguide. In order that the experi¬ 
mental procedure yields satisfactory results, attenuation 
should be taken into account and controlled. Therefore, 
to minimize the propagation losses, the waveguide has 
to be relatively short, of the order of one meter. Ra¬ 
diation losses, on the other hand, are enhanced as the 
number of the scatterers (holes) increases, due to the 
acoustic radiation to the outer environment. An opti¬ 
mal choice which satisfies these experimental objectives 
is schematically illustrated in Fig. [^ and described be¬ 
low. Nevertheless, the emerging (acoustic) structure can 
be considered as a weakly lossy system. Let us note 
that the periodic/aperiodic distribution of strong losses 
can lead to interesting phenomena like formation of an- 
tibandgap [441145] , 

The final aperiodic CLS setup is comprised of seven 
scatterers of two different potential strengths Ui = 
7.8886 and U 2 = 12.3414 corresponding to holes of radii 
Ti = 0.009 m and r 2 = 0.012 m, respectively (for the 
calculation of potential strengths, we use Eq. @). Our 
choice to use a limited number of scatterers in order to 
minimize losses, results in the identification of a single 
decomposition in reflection symmetric parts (defining the 
domains I?i, 7 ) 2 ), which are depicted by the solid arcs in 
Fig.^d). Accordingly, the correspondence between dif¬ 
ferent decompositions in reflection symmetric parts and 
s-PTRs dictates that from this setup only one s-PTR can 
emerge. Selecting the s-PTR frequency to be fs = 450 Hz 
and implementing the corresponding construction prin¬ 
ciple, we find that the distances between the scatterers 
are Li = 0.07 m, L 2 = 0.2671 m and L 4 = 0.1925 m. 

Note that the s-PTR is not affected by the variation 
of the length L 3 separating the two reflection symmetric 
parts of the decomposition. This free parameter will be 
utilized for the construction of an additional a-PTR in 
the following manner: By identifying two parts of the 
device i.e. TT-i, 77 . 2 , we calculate independently their 
transmission spectra and seek for frequencies with com¬ 
mon transmission coefficient. For the case studied here 
this condition is fulfilled at fa = 820 Hz. Fixing this 
frequency and varying the distance L 3 we find that the 
a-PTR occurs when L 3 = 0.1047 m. 
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C. Invariants of broken discrete symmetries in 
systems with losses 

In [^, it has been shown that in the case of scatter¬ 
ing systems characterized by a globally broken discrete 
symmetry, a nonvanishing non-local, spatially invariant 
current Q exists that characterizes the symmetry break¬ 
ing. This analysis, however, concerned lossless idealized 
systems. In experimental, realistic setups losses are in¬ 
evitable and their theoretical account is necessary to ac¬ 
curately describe or predict the expected outcome. In the 
following, we develop a formalism which systematically 
describes the breaking of discrete symmetries in systems 
with attenuation and we elucidate how the quantity Q is 
affected. 

Losses may emerge either during the propagation or 
locally at the scatterer points due to acoustic radiation 
to the outer environment and viscothermal effects. Prop¬ 
agation losses for structures of small length can safely be 
neglected when compared with the attenuation emerging 
from the local scatterers. This means that the impedance 
of the hole becomes a complex function and therefore the 
potential obtains the following form: 

U{x) = Ur{x) + iUi{x). (15) 
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FIG. 3. (Color online) (a) Schematic of an infinite sys¬ 
tem with global reflection and translation symmetry, (b) The 
global reflection symmetry is broken due to the asymmetric 
asymptotic conditions. However, the setup is still reflection 
symmetric, (c) The reflection symmetry of the setup is explic¬ 
itly broken under the presence of two defects. Nevertheless, 
it is retained locally in the restricted spatial domains ©i, 1^2 ■ 


The condition that allows the applicability of the the¬ 
ory is that both the real and the imaginary parts of 
the potential U{x) obey the linear symmetry transform 
U{x) = U{F{x)) for F{x) = X = ax + r]. Depending on 
the value of the parameters a and 77 , F{x) describes a re¬ 
flection at a point a (a = —1, rj = 2a) or a translation by 
a length f} [a = \, j] = f3). The following analysis holds 
for both symmetry transformations whereas in this work 
we focus exclusively on reflection symmetric potentials, 

Ur{x) = Ur{x), 

Ui{x) = Ui{x). (16) 

The potential in Eq. © which is used in this work satis¬ 
fies both conditions. 

To obtain Q, we multiply Eq. ([^ withp(a;) and its par¬ 
ity transform with p{x). Then we subtract the resulting 
equations and due to the symmetry expressed in Eq. (16) 


the potential terms cancel and the resulting expression 
can be written as a total derivative. This, in turn, leads 
to the intriguing result that Q remains spatially invariant 
even in the presence of attenuation, 

Q= Yi = const, (17) 

indicating the strength of the constrain ts w hich are im¬ 
posed by the symmetry. Dividing Eq. 0 by p^{x) we 
obtain, 


Q 1 f ap{x)p'(x) - p{x)p'(x) 


p‘^{x) 2i 


p‘^{x) 


or 


Q 1 (p{x)\' 


p^{x) 2i \p{x) 


(18) 


(19) 


By integration we find 

p{x) = cp{x) - 2iQ p{x) / (20) 

Jx P iO 


where c = p{C)/p{C)- Equation (20) constitutes a map¬ 
ping for the pressure field, linking p{x) from the source 
part to the target part of the symmetry domain, which 
are related via the symmetry transformation x ^ x. Ob¬ 
viously, this transform maps [a:, C] to [a;,C]- This prop¬ 
erty is of fundamental significance, since it generalizes 
the well-known Bloch and parity theorems, not only for 
systems with globally broken reflection symmetry [26) . 
but also when they sustain losses. The latter renders the 
constancy of Q directly observable in realistic and exper¬ 
imentally realizable lossy wave scattering setups. The 
experimental verification of this result will be presented 
in the following section, focusing on the case of reflection 
symmetry. 

The importance of the invariant current Q can be 
perceived within the systematic description of symme¬ 
try breaking. In a globally symmetric system the cor¬ 
responding symmetry is satisfied V a; G K. A reflection 
(or translation) symmetric system extending infinitely, 
as the one illustrated in Fig. (a) belongs to this class. 
A vanishing Q value indicates this case. In scattering 
systems, however, the global symmetries appear broken 
and Q is the key tool to track the mechanism of the sym¬ 
metry breaking. Here it is possible that the structure 
possesses a global symmetry, whereas this symmetry is 
globally broken due to the imposed asymptotic condi¬ 
tions. In this case, which is illustrated in Fig. (b), the 
invariance of Q is retained globally, in the same sense 
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that the energy current is constant in the lossless case. 
For the typical scattering asymptotic conditions, 

Pl(x) = x - Xq < 0 

PB.{x) x-XL>0 (21) 


with pl and pn being the pressure field on the left and 
right side of the setup, respectively, L = xl — xq being 
the length of the setup (see Fig. (b)) and r and t corre¬ 
sponding to the pressure field reflection and transmission 
amplitudes, it is straightforward to show that 


IQI = kT^, for reflection symmetry (22) 

\Q\ = k{RT)^ , for translation symmetry (23) 

with T = |tp and R = |rp being the pressure field trans¬ 
mission and reflection coefficients for the lossy system, 
respectively. Particularly for the reflection symmetry 
transformation, Eq. (22) has a similar form as the en¬ 
ergy current J = fcTiossiess for the lossless analogue, dif¬ 
fering only in the power of the transmission coefficient. 
It is worth mentioning that despite the fact that in a 
system with attenuation the energy current J ceases to 
be constant, the symmetry induces a constant current 
regardless of the losses. 

A more complete breaking of the underlying symmetry 
however takes place if we do not only introduce an asym¬ 
metry of the scattering asymptotic conditions. In this 
case the potential itself globally does not obey the corre¬ 
sponding symmetry transformation but the symmetry is 
present locally i.e., in some spatially limited domains. An 
indicative such case is shown in Fig.|^(c) where the reflec¬ 
tion symmetry is globally broken due to the presence of 
defects, nonetheless it is locally retained in the domains 
I?i, I? 2 - Whereas here there is no spatially invariant 
current Q for the total extent of the device, there exist 
two different currents Qi, Q 2 spatially invariant within 
the symmetry domains I?i, I? 2 , respectively. Breaking 
further the symmetry would lead to smaller regions of 
constant Qs. In the extreme case of a complete break¬ 
ing of the corresponding symmetry i.e. strongly disor¬ 
dered systems without any persisting local symmetry, Q 
ceases being constant. In this sense, the landscape of Q 
provides a systematic tracking of the breaking and the 
corresponding remnants of a symmetry. 


IV. RESULTS 
A. Transmission spectrum 

Having determined the geometric characteristics of the 
waveguide, we study experimentally the local symmetry 
induced characteristics at the PTR frequencies. Figure]^ 
shows the comparison between the theoretical lossless, 
the theoretical lossy and the experimental transmission 
spectra. The theoretical transmission, both for the loss¬ 
less and lossy case, are obtained by the transmission ma¬ 
trix method (see Appendix for details). The experi¬ 
mental transmission spectrum is obtained by calculating 


the transmission matrix elements following the procedure 
described in |33j and using an anechoic termination. In 
addition, due to the asymmetry of the structure, two dif¬ 
ferent measurements (source at the right of the structure, 
anechoic end at the left for the first measurement and the 
reverse position for the second) are performed. 

The experimental transmission spectrum is in very 
good agreement with the lossy theory. Small discrep¬ 
ancies between theory and experiment can be attributed 
to the reflections at the termination end, to a non-perfect 
calibration of the microphone and to unavoidable leaks 
of the system. The obvious effect of losses is the de¬ 
crease of the PTR peaks at a percentage of 40 — 50%. 
Note also that in the lossy case the peaks corresponding 
to PTRs are slightly shifted and appear at frequencies 
fs = 446 Hz and fa = 815 Hz instead of the theoreti¬ 
cally predicted values of fs = 450 Hz and fa = 820 Hz, 
respectively. 


B. Pressure profiles 


Figures 0 (b)-(c) illustrate the magnitude of the pres¬ 
sure field, ^(a;) I, and its corresponding phase 9p for fs, 
fa and fn = 659 Hz as indicated in Fig. (a) by the 
dashed lines, where /„ corresponds to a non-PTR state. 
Experimental measurements of the magnitude and the 
phase of the pressure field at these frequencies were per¬ 
formed by a 1/4 inch microphone (GRAS Type 26AC) 
placed within the tube and positioned at lattice points 
oi dx = 1 cm, see Eig. [^a). 

The theoretically anticipated results (with losses in¬ 
corporated) are obtained by the transmission matrix (see 
Appendix]^ for detailed calculations) using anechoic ter¬ 
mination conditions, namely an impedance termination 
equal to the acoustic characteristic impedance of the tube 
as boundary condition. Theoretical and experimental re¬ 
sults, both for \p{x)\ and for 0 p are in very good agree¬ 
ment. We also observe that the magnitude of the pressure 
field at the s-PTR shown in Fig. Bb) follows the symme¬ 
tries of the mirror symmetric domains which constitute 
the corresponding decomposition, as expected from the 
lossless theory m- The small deviations from the exact 
symmetry are attributed to the losses sustained by the 
system and to the small back propagation effects due to 
imperfections of the anechoic end. 

The spatial profiles of the magnitude of the pressure 
field indicate an additional interesting property of the 
CLS class of devices, which concerns the localization of 
the field. Particularly, Fig. (b) indicates that at the 
s-PTR frequency the localization of the pressure field is 
enhanced in the first mirror symmetric domain compared 
to the second. Similar results have been theoretically 
observed in mills, evincing the potential usefulness of 
this class of structures for the control of localization. 
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FIG. 4. (Color online) (a) Transmission spectrum of the ape¬ 
riodic waveguide: lossless theory (magenta short dash curve), 
theory with losses (purple long dash curve) and experiment 
(green solid curve). The vertical (solid) lines indicate the 
frequencies we focus on: s-PTR frequency (/s), a-PTR fre¬ 
quency {fa) and non-PTR (/n). (b) Pressure field magni¬ 

tude \p{x)\ across the waveguide at the frequencies fs,fa,fn 
marked in (a). The solid and dashed vertical lines show the 
positions of the centers of the holes with two different radii 
and thus different scattering strength. The shaded areas de¬ 
pict the two domains of reflection symmetry (Di, D 2 ) to 
which the setup can be completely decomposed, (c) Phase 
of the pressure field 9p across the waveguide for the three 
selected frequencies /«, /„ and fa. 


C. Observation of the symmetry induced invariant 
currents 


Let us experimentally verify the presence of the in¬ 
variant current Q. For reasons of consistency we choose 
for the analysis the same selected frequencies, fs, fa and 
fn to calculate the current Q. We follow two different 
procedures for this demonstration. 


1. Plane wave fitting method 


This method assumes an one-dimensional sound prop¬ 
agation. Using the experimental pressure field data, we 
perform a fitting procedure with a plane wave function 
of complex wavenumber k = kn — ikj given by the ex¬ 
pression 

f{x) = ~ ^R^)) 

+ \B\e'"‘^ cos{y - {Sb - kRx)). (24) 

We fit simultaneously the real and imaginary parts of 
the experimental pressure field using |yl|, \B\, Oa, 0b, ku 
and fc/ as shared parameters for the two datasets (real, 
imaginary part). When the parameter e is zero, the fit¬ 
ting function corresponds to the real part of the field, 
whereas for e = 1 Eq. (24) corresponds to the imaginary 


part. Then, keeping the complex wave number fixed, we 
find the plane wave parameters A and B. 

However, around the holes, the scattered waves have 
an evanescent part with 3D characteristics which locally 
alter the plane wave description. Thus, we ignore the ex¬ 
perimental points residing in the areas which extend one 
millimetre on either side of each hole and where a devi¬ 
ation from the plane wave description occurs. The re¬ 
maining experimental points comprise five regions -three 
in the domain and two in the domain I? 2 - where the 
plane wave approximation is accurate. Thus, we perform 
a separate fitting procedure in each of these five regions. 
The invariant currents Q are determined by using the 
results of the fitting procedure in Eq. Figures [5] 

(a), (b) illustrate \Q\ and the respective phase 9q for the 
three selected frequencies. 


2. Mapping method 


The alternative way to calculate Q uses directly the ex¬ 
perimental data, without excluding any regions or using 
fitting procedures. Having measured the pressure field in 
the total extent of the acoustic structure, we can directly 
compute Q via the mapping relation. 


Q 


p{x) — p{2a — x) 
2 * 


(25) 


which results from Eq. ( [^ . Figures (c), (d) illus¬ 
trate the corresponding results for \Q\ and 0q, as emerged 
from the latter. The integral in Eq. ( [20| becomes zero at 
the mirror axis position of each domain and the expres¬ 
sion becomes indeterminate. A more accurate approach 
of this vicinity would require more experimental data 
points. Therefore we have ignored a small area around 
the domain center. This also causes the observed increas¬ 
ing trend as the center of the domains is approached. 
The small fluctuations of the values of \Q\ and 9q are 
inevitable since the experimental points deviate from the 
exact plane wave which would be the solution of Eq. (|^ . 
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FIG. 5. (Color online) Schematic of the setup, illustrating the spatial profile of \Q\ and 6q. The solid and dashed lines 
correspond to the positions and strengths Ui and U 2 of the point scatterers, respectively. The coloured areas delimit the 
domains Di , ID 2 from the start of the first hole to the end of the final hole. The magnitude and phase of the invariant current 
Q are calculated for the selected s-PTR (/^ = 446 Hz), a-PTR {fa = 815 Hz) and non-PTR (/„ = 659 Hz) frequencies, (a) and 
(b) show the magnitude and the phase of Q calculated with the method of fitting the experimental data in each potential free 
area with the proper plane wave. Accordingly, (c) and (d) show \Q\ and 6q calculated using the mapping relation (Eq. (20l). 


3. Remarks 


Comparing the results of the two approaches it is ob¬ 
vious that they are in very good agreement and the con¬ 
stancy of Q in the symmetry domains becomes evident, 
verifying the theoretical prediction. In particular, since 
our structure is composed of two local symmetry sub- 
domains Vi, I? 2 , according to the lossless theory [30] 
the magnitude of the invariant currents \Q\ in each do¬ 
main at the PTR frequencies should be the same. At any 
other non-PTR frequency, the magnitudes of the invari¬ 
ant currents Q is constant in each domain of the selected 
decomposition, but they differ from each other. 

In the first place, both for theory and experiment (see 
Figs. (a) and (c)), we observe that in the presence 
of attenuation, |Q| is different at the two decomposition 
domains even at the PTR frequencies, thereby deviating 
from the lossless theory. Using the theoretical transmis¬ 
sion matrix method and tuning the total amount of losses 
from zero up to the expected value, we observe that as at¬ 
tenuation increases, the larger the difference between \Q\ 
in Vi and I? 2 , becomes. On the other hand, for a given 
amount of losses which is kept hxed, by changing the fre¬ 
quency, we observe that this difference is minimized at 
the s-PTR frequency, being in accordance to the lossless 
theory |5S]. 

An additional feature which is observed in the plane 
wave fitting method and is worth-mentioning, concerns 


the small gaps which appear in the experimental data 
within the domain Vi (resp. ID 2 ) and particularly be¬ 
tween the hrst (or third for Vi) region and the second. 
This can be explained by the three-dimensional nature 
of the problem in the vicinity of the holes. This feature 
alters the solution from one region to the other, leading 
to disconnected plane wave solutions for each one of the 
three regions, instead of the exact, continuous plane wave 
for the whole first domain, as in the case of the theoreti¬ 
cal results. This difference in the plane wave solutions is 
therefore translated into a gap for \Q\ and 0 q. 

Concluding this subsection, we point out that a com¬ 
plete constancy of Q would be feasible if the pressure field 
in every symmetry domain was a perfect plane wave, sat¬ 
isfying the ID Helmholtz equation. The deviation of the 
experimental curve from the exact functional form of the 
plane wave results in a gap in domain when the fitting 
procedure is followed or in small fluctuations when the 
calculation is performed via Eq. (201. 


D. Field mapping between symmetry related areas 

In Ref. |55|, it was demonstrated how the spatially in¬ 
variant current Q constitutes the key tool for tracking 
the symmetry breaking procedure in lossless, aperiodic 
systems. Here, having already identified the constancy 
of Q in the presence of losses, we experimentally ver- 
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X (m) 



FIG. 6. (Color online) Real and imaginary parts of the pressure field p{x) for the selected s-PTR {fa = 446 Hz), a-PTR 
{fa = 815 Hz) and non-PTR (/„ = 659 Hz) frequencies. The black (large) circles correspond to the experimental field’s spatial 
profile. The green small circles (always coinciding with the black larger circles) also correspond to the experimental data and 
are used to indicate the experimental data points (source domain in white) which are mapped through Eq. (201 to the magenta 
diamonds (target domain in magenta). 


ify the generalization of the parity theorem in the case 
of complete local reflection symmetry in the presence of 
losses. To this aim, we utilize the aperiodic setup con¬ 
sidered above (see Fig. and use the measured Qs and 
pressure fields for the frequencies fs, fa and /„. The 
mapping between the symmetry related regions is pro¬ 
vided by Eq. (1^. 

In Fig. [^the real and imaginary parts of the pressure 
field are illustrated for the selected s-PTR (/^ = 446 
Hz), a-PTR {fa = 815 Hz) and non-PTR (/„ = 659 
Hz) frequencies. To investigate the validity of the above 
mapping relation expressed in Eq. (20), we employ as 
input parameters the experimental pressure field at the 
left side of each mirror symmetry axis, represented by 
the green (small) circles and the Q values obtained by 
the fitting procedure, presented in the previous section 
(Fig. [5] (a), (b)). Note that we use the Q values obtained 
by the fitting method in order to have an independent 
mapping result. The resulting mapped points, denoted 
by the magenta diamonds, are in very good agreement 
with the experimental data. 


V. CONCLUSIONS 

This is the first experimental implementation and ver¬ 
ification of the concept of local symmetries and the con¬ 
struction principle proposed in Refs. [291 ED] for aperi¬ 
odic acoustic waveguides with prescribed perfect trans¬ 
mission resonances. The experimental observation re¬ 


veals the universal character of this concept, which has 
been theoretically predicted to occur in quantum me¬ 
chanical aperiodic structures |29j and photonic multilay¬ 
ered structures [30], but is now experimentally verified 
in a completely different, namely acoustic, system. In 
particular, we extend the theory by incorporating losses 
and ascertain that the symmetry induced invariant cur¬ 
rent Q preserves its spatial constancy in the presence of 
attenuation. The fact that broken discrete symmetries, 
retained at spatially limited domains, induce invariant 
currents, has been experimentally verified here for the 
case of locally reflection symmetric setups. In turn, we 
showed how Q generalizes the Bloch and parity theo¬ 
rems for scattering systems with broken global transla¬ 
tion and reflection symmetry, respectively, by defining 
a mapping of the pressure field between the symmetry 
related regions. Subsequently, we experimentally verified 
the validity of this generalization for the case of reflection 
symmetries. It has been demonstrated how Q constitutes 
the key-tool towards a systematic description of discrete 
symmetry breaking. In this context, we showed that in a 
globally reflection symmetric setup which sustains losses, 
IQ I provides the lossy analogue to the energy current J of 
the lossless case. Finally, the construction procedure for 
the design of acoustic waveguides with prescribed perfect 
transmission properties and the subsequent experimental 
implementation reveals how the class of CLS materials 
could pave the way towards the design of materials with 
very specialized filtering capabilities. 
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Appendix A: Length corrections due to the hole where 


For fcr„ ^ 1, Refs. m and [35] give the following 
expressions for, respectively, the length correction due 
to radiation inside the principal waveguide and due to 
radiation to outer environment, respectively: 


T = 


'Nh-1 


n 


n—1 


M 


(h) 

^ \ U, 


with 


= 0.82(1-0.235e„-1.32e2+1.54e3-0.86e^)r„ (Al) 


= 


COs{kdn) 
—ijZc. sin(A:dri 


—iZc sin(A;d„ 
cos{kdn) 


and 


and 


= (0.82- 0.47(°«)e„)r„, (A2) 



1 0\ 
1/Zn 1 ) • 


(B2) 


(B3) 


(B4) 


with e„ = r„/i?. 


Appendix B: Transmission Matrix Method 


Applying the transmission matrix method, the propa¬ 
gation through a finite system made of Nh holes can be 
expressed as follow |33| 


{ Pzn 


= T 


Pout 

Uout 


(Bl) 


and describe the propagation through a 

waveguide of length dn = Xn+i — Xn and through a 
branched scatterer with input impedance respec¬ 
tively. Pin iJJin) and Pout iUout) represent the pressure 
(volume velocity) at the input and output of the system, 
respectively. This transmission matrix method allows to 
determine directly the pressure field inside the structure 
and the transmission coefficient t of the aperiodic system 

by [33] 


2 

Til + TvilZc. -b Ti\Zc -b T22 


(B5) 
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